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l  introduction 

Scattering  phenomena  occur  universally  for  short  period  (frequency  ^  1  Hz)  seismic 
waves  observed  at  the  earth's  surface.  One  of  the, major  manifestations  of  scattering  Is 
the  presence  of  coda,  seismic' energy  following  a  major  P  or  S  arrival  that  cannot  be 
explained  by  propagation  in  a  layered  structure.  In  this  project  coda  is  being  examined  by 
two  methods.  The  first  method  is  to  use  the  NORESS  array  in  southern  Norway  to  examine 
teleseismic,  regional  and  local  coda.  In  preliminary  investigations  carried  out  last  year  on 
teleselsmlc  P  code  (Dainty,  1986),  a  promising  analysis  method  ("wavenumber  spectrum") 
was  Identified.  In  this  report  we  discuss  some  tests  that  have  been  made  to  characterize 
certain  biases  and  limitations  of  the  wavenumber  spectra  noted  in  the  preliminary  investi¬ 
gation.  The  second  method  we  are  using  Is  ultrasonic  modeling;  preliminary  work  on  building 
models  to  produce  coda  in  the  laboratory  Is  presented  together  with  a  discussion  of  some 
appropriate  investigations  that  might  be  carried  out. 

The  objective  of  the  study  Is  to  understand  the  coda  so  that  it  may  be  used  to  esti¬ 
mate  yield  (Bullitt  and  Cormier,  1984;  Gupta  et  a i„  1986),  measure  attenuation  (Aki  and 
Chouet,  1976)  and  as  a  possible  discriminant  between  earthquakes  and  explosions 
(Dainty,  1986).  At  present  these  objectives  suffer  from  a  lack  of  understanding  of  the 
components  that  make  up  coda  and  uncertainty  concerning  the  limits  of  validity  of  the  sin¬ 
gle  scattering  theory  used  to  analyse  coda. 
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D.  WAVENUMBER  SPECTRA 


{Summary  of  Previous  Work 

The  data  taken  at  NORESS  may  be  considered  as  a  sampling  in  time  and  two  space 
dimensions  of  a  wavefleld.  Since  harmonic  waves  are  fundamental  solutions  of  the  wave 
equation,  it  is  often  useful  to  perform  a  three  dimensions!  spectral  analysis  to  find  power 
as  a  function  of  frequency  /  and  two  components  of  (linear)  wavenumber,  ssy 
wavenumber  North  (kg)  and  wavenumber  East  (kg).  A  contour  plot  of  such  a  transform  at  a 
frequency  of  3.6  Hz  against  kg  and  kg  is  shown  in  Figure  1  for  a  teleseismic  event  shown, 
along  with  the  coda  analysis  window  used,  In  Figure  2.  From  such  a  plot  we  may  begin  to 
break  down  the  coda  into  constituent  parts  corresponding  to  waves  scattered  in  different 
places.  Two  well  known  criteria  may  be  used.  One  uses  the  azimuth  of  a  wave  measured 
clockwise  from  North,  given  by 

tarn?  r  kE/ kN  (1) 

The  second  criterion  is  to  use  the  apparent  velocity  Va  given  by 

Va  -  f/k  (2) 

where  the  total  horizontal  wavenumber  k  is  given  by 

k2  =  kg2  +  kg 2  (3) 

The  apparent  velocity  *8  related  to  the  velocity  of  the  surface  material  V,  and  the  angle 
of  incidence  of  the  wave  to  the  vertical  i  by 

Va  -  7,/slni  (4) 

In  the  coda  i  is  not  In  general  known,  and  another  relation  from  (4)  must  be  used  Instead, 
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MONOCHROMATIC  F-K’ANALYSIS 

CONVENTIONAL  RESOLUTION 

LINEAR  DEPENDENCY  BETWEEN  CONTOURS 


PREFILTER  FREQUENCY  INTERVAL:  0.0  -  0.0  HZ 

COOA 


ANALYSIS  r«.-.-£NCY:  5.6  HZ 


figure  i .  Frequency-wavenumber  spectrum  for  event  and  window  of  figure  2.  Analysis 
frequency  Is  3.2  Hz,  KX  is  wavenumber  East,  KY  wavenumber  North,  units  are  1  /km.  From 
Dainty  (1986). 
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For  Investigations  of  coda,  the  Interesting  ranjj  i  of  ?„  Is  the  crustal  range  3  -  9  km/s;  any 
component  with  an  apparent  velocity  faster  4jn.n  9  km/s  may  be  classified  as  teleseismlc. 

From  (1),  (2)  and  (3)  lines  passing  outwards  from  the  origin  on  Figure  1  are  linos  of 
constant  aJdmuth  while  circles  cantered  on  the  origin  are  lines  of  constant  total 
wavenumber  and  also  lines  of  constant  apparent  velocity.  A  line  of  constant  azimuth  and 
two  circles  of  constant  phase  velocity  are  shown  on  Figure  1 ;  teleseismic  P  waves  will 
have  faster  apparent  velocities  than  8  Km/s  while  Lg  waves  will  have  velocities  close  to 
3.6  km/s.  In  coda  studies  it  Is  often  useful  to  know  the  total  power  as  a  function  of 
azimuth  or  total  wavsnumber  (l.e.,  apparent  velocity).  To  obtain  this  Information  the 
frequency-wavenumber  spectrum  shown  in  Figure  1  must  be  Integrated  along  a  line  of  con¬ 
stant  azimuth  or  constant  total  wavenumber.  Appropriate  formulations  for  doing  this  will  be 
given  In  a  forthcoming  publication.  Figure  3  shows  the  resuit  for  azimuth  and  Figure  4  for 
wavenumber.  So  far,  the  wavenumber  spectrum  as  shown  In  Figure  4  has  been  used  in 
coda  analysis. 

Figure  4  Illustrates  some  of  the  strengths  and  also  one  of  the  limitations  of 
wavenumber  spectra.  Note  that  there  are  two  peaks  superimposed  on  a  linear  trend.  The 
two  peaks  Indicate  that  there  is  a  preponderance  of  energy  at  two  phase  velocities  In  the 
signal,  a  high  phase  velocity  of  1 8  km/s  and  a  low  phase  velocity  of  4.3  km/s.  The  high 
phase  velocity  peak  may  also  be  seen  In  Figure  1  (the  peak  Inside  the  8  km/s  circle) 
where  it  Is  clear  that  the  energy  Is  coming  In  with  a  definite  azimuth  (~  75°)  that  is  in  fact 
the  source  azimuth,  Indicating  near  source  scattering.  The  high  phase  velocity  indicates 
that  the  energy  has  travelled  from  the  source  as  P.  A  comparison  of  Figure  4  and  Figure  1 
Indicates  that  the  low  velocity  peak  is  made  up  of  the  sum  of  several  peaks  with  differing 
azimuths,  nene  of  them  the  source  azimuth.  This  diffuse  component  of  the  coda  Is  Inter¬ 
preted  as  being  locally  scattered  within  the  crust  near  the  receiver,  since  it  has  the 
apparent  velocity  of  S  waves  Incident  needy  horizontally.  Finding  such  "diffuse"  com- 
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ponents  in  the  coda  has  been  a  major  use  of  the  wavenumber  spectrum. 


One  of  the  difficulties  encountered  in  the  use  of  wavenumber  spectra  has  been  the 
presence  of  linear  trends  such  as  that  labelled  "white  noise"  In  Figure  4.  The  term  is  used 
because  a  constant  ("white  noise")  trsquency-wavenumber  spectrum  will  produce  such  a 
trend  due  to  the  increasing  circumference  of  the  circle  of  integration  as  the  totai 
wavenumber  increases.  This  problem  increases  with  increasing  frequency  if  the  full  array 
is  used.  At  other  problem  shows  up  at  low  frequencies  as  illustrated  by  Figures  5  and  6, 
which  present  analyses  of  the  event  shown  in  Figure  2  but  at  1 .6  -  1 .8  Hz  instead  of  3.6 
Hz.  Figure  5  shows  the  frequency-wavenumber  spectra  for  a  window  around  first  P  and  for 
the  coda  20  -  100  sec  after  first  P.  Both  show  a  prominent  peak  with  the  source  azimuth 
and  a  high  apparent  velocity  -  15.6  km/s  for  first  P  end  11.7  km/s  for  the  coda.  The 
wavenumber  spectrum  for  the  same  two  windows  is  shown  in  Figure  6.  The  first  P  window 
shows  a  single  major  peak  but  the  peak  is  at  a  lower  apparent  velocity  (9  krn/s)  than  the 
peak  In  Figure  5,  i.e.,  there  is  a  bias  in  the  apparent  velocity  determined  from  the 
wavenumber  spectrum.  The  coda  wavenumber  spectrum  has  been  interpreted  as  consist¬ 
ing  of  two  peaks  (Dainty,  1985)  as  Indicated  in  Figure  6.  The  peak  with  the  higher 
apparent  velocity  Is  Interpreted  as  corresponding  to  the  peak  in  Figure  5,  but  again  there 
is  a  bias  towards  a  lower  apparent  velocity  (7.2  km/s).  The  other  peak  corresponds  to  the 
diffuse  local  component  of  the  coda  and  is  not  evident  In  Figure  5,  presumably  because  of 
the  spread  in  azimuth. 

1 

K  Coherence  and  Bias  in  Wavenumber  Spectra 


The  two  problems  noted  above  may  be  explained  by  the  known  characteristics  of  the 
array  and  seismic  signals  received  at  It.  In  this  section  we  present  soms  work  on  the  use 
of  subarrays  to  minimise  the  "white  noise"  linear  trends  observed  at  high  frequencies,  and 
work  on  the  use  of  synthetic  signals  to  quantify  the  apparent  velocity  bias  observed  at 
low  frequencies. 


TT 


Wavenumber  kg 


Figure  6.  Frequency-wavenumber  spectra  of  the  event  shown  In  Figure  2  at  a  frequency 
of  1.6  Hz.  (Top)  Five  second  window  around  first  P.  (Bottom)  Wndow  between  20  and 
100  seconds  after  first  P  (coda),  calculated  as  average  of  non-overlapping  6  second  win¬ 
dows.  From  Dainty  (1986). 
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Figure  6.  Wavenumber  spectra  of  the  event  shown  In  Figure  2  at  a  f;«oqu*ncy  of  1.8  Hz 
Wndows  as  for  Figure  5.  From  Dainty  (1986). 


Obe> possible  explanation  for  the  white  noise  phenomenon  is  a  lack  of  correlation  of 
the  signal  at  different  seismometers  of  the  array.  Uncorrelated  signals  would  appear  as 
white  noise  on  a  wavenumber  spectrum  because  they  cannot  be  expressed  as  travelling 
waves.  Such  a  lack  of  correlation  at  distances  comparable  to  the  inter-seismometer  spac¬ 
ing  has  been  observed  in  studies  of  noise  at  NORESS.  Figure  7  from  Bungum  et  a l.  (1985) 
shows  noise  correlation  against  distance  for  the  NORESS  site  for  four  frequency  bands.  If 
the  noise  can  be  considered  to  be  a  wavefield  consisting  of  a  mixture  of  travelling  waves, 
the  results  will  also  apply  to  seismic  signals,  especially  in  the  coda.  From  Figure  7,  at  a 
fixed  seismometer  separation,  500  m  say,  the  correlation  declines  with  increasing  fre¬ 
quency.  A  map  of  the  NORESS  array  is  shown  In  Figure  8,  and  It  may  be  seen  that  the 
outermost  ring  of  seismometers  has  a  minimum  separation  of  about  500  m  from  other  array 
locations.  From  Figure  7  a  deterioration  of  the  wavenumber  spectrum  at  frequencies 
greater  than  2  Hz  would  be  expected  if  the  outer  ring  Is  used;  this  Is  observed.  The  effect 
should  be  more  pronounced  at  higher  frequencies;  this  also  is  observed. 

If  this  is  in  fact  the  cause  of  the  white  noise  problem,  It  should  be  possible  to  mitigate 
It  by  eliminating  from  the  calculation  of  the  wavenumber  spectrum  array  elements  that  are 
far  from  all  other  array  elements.  Tests  have  been  made  of  subarrays  from  which  the  outer 
ring  has  been  eliminated  and  from  which  the  two  outer  rings  have  been  eliminated  (Figure 
8).  Figure  9  shows  the  seismogram  at  the  center  seismometer  for  a  Iccal  quarry  blast 
south  of  NORESS  with  an  analysis  window  around  first  P  indicated.  Figures  10  -  15  are 
frequency-wavenumber  and  wavenumber  spectra  for,  in  order,  the  full  array  (taken  at  7.2 
Hz,  close  to  the  peak  signal  frequency),  the  array  minus  the  outer  ring  (subarray  A),  and 
the  array  minus  the  outer  two  rings  (subarray  B).  While  the  frequency-wavenumber  spec¬ 
trum  using  the  full  array  (Figure  1 0)  does  pick  out  the  appropriate  azimuth  and  phase  velo¬ 
city  of  the  arrival,  the  wavenumber  spectrum  (Figure  11)  is  heavily  contaminated  with 
noise,  although  there  is  a  peak  at  the  correct  wavenumber  (1  km"1)  The  results  for  subar¬ 
ray  A  (Figures  12  and  13)  and  subarray  B  (Figures  14  and  15)  have  considerably  less 
noise,  but  besides  a  peak  at  the  expected  wavenumber  of  1  km’1  subarray  A  shows  a  peak 
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Figure  7.  Observed  correlation  of  noise  between  array  sensors  as  a  function  of  distance 
between  them  taken  using  a  temporary  NORESS  array  of  21  seismometers.  Data  was  col¬ 
lected  on  August  1 2,  1 083.  Frequency  band  analysed  Is  shown  on  each  plot;  fell  plots  use 
the  same  data.  From  Bungum  ef  ol.  (1985), 
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Rgure  1 0.  Frequency-wavenumber  spectrum  of  event  and  window  of  Figure  9  at  a  fre¬ 
quency  of  7.2  Hz  using  the  fun  array. 
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Rjjura  1 2.  Frequency-wevenumber  spectrum  of  event  and  window  of  Figure  0  *t  »  fre¬ 
quency  of  7.2  Hz  using  subarray  A. 
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Figure  1 5.  Wavenumber  spectrum  of  even' 
Hz  using  subarray  B. 


at  a 'Wavenumber  of  2.4km"1(apparent  velocity  3  km/s;  Figure  13)  that  Is  not  evident  on 
plflure'15;(subarray  B); 

The  peak  at  2.4  km"1  appears  to  be  an  artifact  caused  by  the  sidelobe  pattern  of 
subarray  A;  a  similar  peak  is  seen  on  Figure  T1  using  the  full  array.  Figures  16-18  show 
the  beam  patterns  for,  in  order,  the  full  array,  subarray  A,  and  subarray  8.  Figures  16  and 
17  (the  full  array  arid  subarray  A)  have  the  first  sldelobes  at  wavenumbers  of  1.9  km”1 
from  the  central  peak,  while  Figure  1 8  (subarray  B)  has  the  first  sidelobes  at  3  km-1. 
These  differences  reflect  the  smaller  aperture  and  the  smaller  average  distance  between 
seismometers  of  subarray  B.  Because  the  sldelobfjs  are  closer  to  the  main  peak  for  the  full 
array  and  subarray  A,  spurious  peaks  due  to  the  sldelobes  will  appear  at  lower 
wavenumbers  than  for  subarray  B.  Tills  is  a  potentially  serious  difficulty  when  searching 
for  the  diffuse  coda  components  mentioned  previously.  In  assessing  the  severity  of  the 
problem,  note  that  the  lowest  apparent  velocity  of  Interest  is  3  km/s  and  that  the 
wavenumber  corresponding  to  this  may  be  calculated  from  (2)  at  any  given  frequency,  if 
the  full  array  or  subarray  A  are  to  be  used  this  maximum  wavenumber  should  be  s  2  km-1, 
while  for  subarray  B  tho  maximum  wavenumber  should  be  <  3  km";. 

The  problem  of  bias  In  the  velocity  determination  using  wavenumber  spectra  hes  been 
examined  through  the  use  of  synthetic  data.  Figure  1 9  shows  a  synthetic  signal  for  the 
center  seismometer  of  the  array.  A  similar  signal  was  placed  on  the  other  seismometers  of 
the  (full)  array  with  the  appropriate  time  shifts  for  a  P  wave  coming  from  due  North  with  an 
angle  of  incidence  of  15°  to  the  vertical.  A  surface  velocity  V,  of  6  km/s  was  assumed, 
corresponding  to  an  apparent  velocity  Va  of  23.2  km/s,  and  a  frequency  of  2  Hz.  Figure 
20  shows  the  frequency-wavenumber  spectrum;  the  correct  velocity  and  azimuth  are 
reproduced  In  the  position  of  the  peak.  Figure  21  shows  the  wavenumber  spectrum.  The 
peak  occurs  at  a  wavenumber  of  0.18  km"1  corresponding  to  a  velocity  of  11  km/s,  con¬ 
firming  the  bias  suggested  previously.  This  bias  Is  caused  by  the  spreading  of  the  peak  In 
the  frequency-wavenumber  plot  due  to  the  finite  array  aperture  combined  with  the 
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Figure  1 8.  Beam  pattern  of  subaiTay  B. 
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Bf/ura  20.  Frequency-wavenumber  spactrum  of  signal  and  window  of  Figure  19  at  a  fra- 
i  quancy  of  2  Hz  using  full  attay. 
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increasing  circumference  of  the  circle  of  Integration*  which  tends  to  weight  larger 
wavenumbers  disproportionately,  especially  If  the  broadened  peak  overlaps  the  origin  of 
the  frequency-wavenumber  plane. 

The  effect  described  above  has  been  evaluated  by  repeating  the  experiment  shown 
in  Figures  19-21  for  frequencies  of  1  and  2  Hz  and  true  apparent  velocities  from  2.0  to 
23.2  km/s.  Results  are  presented  In  Figure  22,  which  plots  the  measured  velocity  using 
the  wavenumber  spectrum  ageinst  the  true  velocity.  Besides  allowing  us  to  estimate 
corrections  to  the  measured  velocity,  Figure  22  demonstrates  that  there  is  a  "saturation" 
phenomenon  at  high  velocities  -  for  1  Hz,  the  maximum  value  of  the  measured  velocity  is 
6.9  km/s  by  this  method,  no  matter  how  fast  the  true  velocity.  Let  the  horizontal 
wavelength  be  defined  as 


A  =  \/k  -  Va/f  (6) 

This  wavelength  describes  the  length  scale  of  the  wave  as  viewed  by  the  array.  If  the 
bias  is  due  to  the  finite  array  aperture,  wa  would  expect  some  measure  of  the  bias  such  as 
the  ratio  R  of  the  measured  to  the  true  velocity  to  depend  only  on  the  ratio  of  A  to  the 
array  aperture  A.  In  Figure  23  the  results  shown  in  Figure  22  are  replotted  as  a  graph  of 
R  against  X/A.  The  points  for  1  and  2  Hz  both  fall  on  the  seme  curve,  indicating  that  the 
bias  Is  indeed  controlled  by  array  aperture.  The  curve  falls  between  two  asymptotes, 
R-  1  forA/j4-*0  and  R  =  2/ (A  /A)  for  A/  A  large.  The  second  asymptote  expresses 
the  saturation  effect  described  above. 

The  results  of  the  Inv&stigattona  Into  bias  have  a  number  of  practical  implications  for 
work  on  coda.  Corrections  to  observed  velocities  may  be  made  using  the  Information 
presented  In  Figures  22  and  23.  The  most  convenient  method  of  doing  this  is  Figure  24, 
which  plots  R  against  X'/A,  where  A*  is  measured  by  substituting  the  measured  velocity 
for  the  true  velocity  in  (6).  The  high  wavenumber  saturation  asymptote  becomes  X'/A  *  2. 
Another  Important  result  of  the  bias  study  is  the  finding  that  It  is  not  possible  to  measure 
velocities  faster  than  a  certain  limit  for  a  given  combination  of  array  aperture  and 
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frequency,  the  practical  limits  are  A/ A  *  2  or  X'/A  *  1.5,  with  the  velocities  then  deter¬ 
mined  through  (6).  At  1  Hz  with  the  . full  array,  for  example,  this  leads  to  a  (true)  velocity 
limit  of  6 km/s,  meaning  that- it  Is. possible  to  distinguish  between  S  waves  and  P  waves 
but  not  between  crustal  P  and  teleseismic  arrivals.  This  is  probably  the  minimum  require¬ 
ment  for  a  useful  array. 

The  final  step  In  this  investigation  is  the  setting  of  frequency  bands  within  which  dif¬ 
ferent  array  configurations  should  be  used.  The  problems  of  coherence  and  sidelobe 
interference  demand  as  small  an  array  as  possible,  while  minimization  of  velocity  bias 
requires  as  large  an  aperture  as  possible  (the  aperture  of  the  full  array  Is  3  km,  of  subar¬ 
ray  A  1 .6  km,  of  aubarray  B  750  m).  A  minimum  velocity  limit  of  6  km/s,  as  discussed  In  the 
jest  paragraph,  Is  used  to  determine  the  aperture.  Then  the  following  combinations  of  fre¬ 
quency  and  array  configuration  are  proposed:  1  -  2.6  Hz,  the  full  array;  2.6  -  6  Hz,  subar¬ 
ray  A;  5  -  1 0  Hz,  subarray  B.  Wavenumber  spectra  should  not  be  used  for  frequencies  less 
than  1  Hz  at  NORESS;  Note  that  the  Investigations  of  Dainty  (1985)  using  the  full  array  In 
the  frequency  range  1  -  2  Hz  meet  the  recommendetions  given  above. 
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III.  ULTRASONIC  MODELLING 

This  section  gives  a  very  brief  report  of  ultrasonic  Investigations  Into  coda.  At  this 
stage,  a  particular  target  is  the  testing  of  the  sjngle  backscattering  theory  of  the  coda 
advanced  by  Aki  (1969)  and.Akl  and  Chduet  (1976).  This  theory  has  become  the  standard 
model  of  the  coda  and  forms  the  basis  of  the  measurement  of  attenuation  and  magnitude 
using  coda,  an  increasingly  important  endeavour.  In  spite  of  this,  there  has  been  no  good 
test  of  the  single  scattering  assumption,  a  crucial  point  in  the  theory,  Dainty  at  al. 
(1984)  suggest  that  In  several  cases  of  local  earthquake  codas  examined  the  assumption 
of  single  scattering  was  violated.  Accordingly,  It  is  proposed  to  test  the  single  back- 
scatter  theory  of  local  coda  by  constructing  two  dimensional  ultrasonic  models  consisting 
of  1  /1 6"  aluminum  plates  with  holes  or  other  shapes  drilled  Into  them  as  scatterers. 

As  a  first  step  one  such  plate  has  been  constructed  (Figure  25).  The  (Identical) 
transducers  are  small  (~  2  mm  diameter)  piezoelectric  crystals;  orientated  as  shown  In  Fig¬ 
ure  26  they  will  produce  and  record  motions  transverse  to  the  plate,  corresponding  to  flex¬ 
ural  modes.  Flexural  modes  show  pronounced  dispersion  If  the  wavelength  Is  greater  than 
•lx  (6)  times  greater  than  the  plate  thickness  but  travel  at  close  to  the  shear  wave  speed 
for  shorter  wavelengths  (Ewing  et  al,  1 967,  p.  286),  with  motion  predominantly  transverse 
to  the  plate.  For  wavelengths  less  than  six  times  the  plate  thickness,  flexural  modes  may 
be  a  good  two  dimensional  analog  of  SH.  This  has  the  advantage  of  making  interpretation 
of  the  results  simpler  because  SH  does  not  convert  to  P  or  SV  upon  reflection  or  scattering. 
Figure  28  shows  a  seismogram  takon  using  the  setup  In  Figure  26.  The  dominant  frequency 
is  about  500  kHz;  since  tha  velocity  of  shear  waves  In  aluminum  Is  8  mm/^3,  the 
wavelength  is  6  mm,  meeting  the  criterion  for  SH  motion;  dispersion  Is  not  evident.  This 
conclusion  is  supported  by  the  travel  time  of  the  dominant  arrival  In  Figure  26,  30  fjs  after 
origin,  appropriate  for  a  shear  wave  travelling  90  mm  (Figure  26).  Some  P  energy  Is  visible 
before  this  time,  but  It  is  small.  We  have  confirmed  that  there  is  little  motion  parallel  to  the 
plate.  Figure  26  may  be  compared  with  a  local  earthquake  seismogram,  Figure  27  (from 
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Figure  25,  Ultrasonic  model  for  coda.  (Top)  Section  view  showing  orientation  of  sensors. 
(Bottom)  Plan  view. 


Dainty  et  ql„  1 984). 

The  final  consideration  in  constructing  the  model  is  the  matter  of  scaling  between  the 
earth  and  the  model.  Obviously,  the  model  is  much  simpler  than  the  earth,  but  some  scaling 
considerations  still  apply.  Since  the  shear  velocity  of  the  crust  is  about  the  same  as 
aluminum  (3  km/s),  wavelengths  scale  Inversely  as  frequencies;  thus  5  Hz  waves  jn  the 
crust,  a  typical  case  from  Dainty  et  al.  (1984),  will  be  equivalent  to  600  kHz  waves  jn  the 
model  if  km  in  the  earth  correspond  to  cm  in  the  model.  In  the  model  shown  in  Figure  26, 
the  scatterers  (holes)  are  5  mm  (about  one  wavelength)  across  ant.-  correspond  to  objects 
of  about  500  m  linear  dimension  In  the  earth,  although  they  are  almost  certainly  stronger 
scatterers  than  found  in  the  earth.  Another  important  factor  in  the  scaling  is  the  length  of 
the  plate,  since  this- controls  how  farwaves  travel  (60  cm)  before  they  can  contaminate 
the  coda  due  to  end  reflections.  This  Implies  that  the  coda  is  uncontaminated  by  end 
reflections  before  170  /is,  l.e.,  most  of  the  record  shown  In  Figure  26.  This  scales  to 
codas  up  to  1 7  s  after  origin  in  the  earth,  covering  many  of  the  codas  considered  In  Dainty 
etal.  (1984). 

The  second  quantity  that  has  to  be  scaled  Is  the  mean  free  path,  which  Is  the  recipro¬ 
cal  of  the  total  turbidity.  Dainty  et  al,  (1984)  found  mean  free  paths  between  10  and 
100  km;  scaling  this  to  the  model  (10  -  100  cm)  means  that  codas  containing  scattered 
waves  that  have  traversed  0.6  -  6  mean  free  paths  may  be  Investigated.  This  appears  to 
be  the  interesting  range.  To  achieve  this,  the  number,  type  and  shape  of  scatterers  will  be 
varied  as  well  a6  the  ultrasonic  frequency  and  the  length  of  the  plate.  A  particular  object 
of  study  will  be  the  transition  from  single  to  multiple  scattering  and  its  effect  on  parame¬ 
ters  dertfred  from  the  coda. 
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